Thermodynamics of dense hadronic matter in a parity doublet model 
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We study thermodynamics of nuclear matter in a two-flavored parity doublet model within the 
mean fleld approximation. Parameters of the model are chosen to reproduce correctly the properties 
of the nuclear ground state. The model predicts two phase transitions in nuclear matter, a liquid- 
gas phase transition at normal nuclear density and a chiral transition at higher density. At flnite 
temperature the pion decay constant exhibits a considerable reduction at intermediate values of 
chemical potential, which is traced back to the presence of the liquid-gas transition, and approaches 
zero at higher chemical potential associated with the chiral symmetry restoration. A "transition" 
from meson-rich to baryon-rich matter is also discussed. 
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1. INTRODUCTION 

Model studies of dense quark matter have suggested a 
rich phase structure of QCD at temperatures and quark 
chemical potentials of order Aqcd- Our knowledge on 
the phase structure however remains limited and the de- 
scription of strongly interacting matter does not reach 
a consensus yet [l|. In particular, properties of baryons 
near the chiral symmetry restoration are poorly under- 
stood. We believe that the realistic modeling of dense 
baryonic matter must take into account the existence of 
the nuclear matter saturation point, i.e. the bound state 
at baryon density po — 0.16 fm~^, like in Walecka type 
models ^] . Several chiral models with pure hadronic de- 
grees of freedom 0, 01 have been constructed in such a 
way that the nuclear matter has the true ground state. 
Other models consider a nucleon as a dynamical bound- 
state of a diquark and a quark Q . 

In chiral models it is mostly assumed that a con- 
stituent quark transforms under chiral rotations as a bare 
quark. According to this naive assignment, dynamical 
chiral symmetry breaking entirely generates a baryon 
mass which thus vanishes at the restoration. Accord- 
ing to an alternative, mirror assignment [a, |7| , dynam- 
ical chiral symmetry breaking generates a mass differ- 
ence between parity partners and the chiral symmetry 
restoration does not necessarily dictate the chiral part- 
ners to be massless. Mirror baryons embedded in linear 
and non-linear chiral Lagrangians have been applied to 
study their phenomenology in vacuum [6|-|8| , nuclear mat- 
ter [^, [i3| and neutron starts [11|, although it remains 
an open question which scenario is preferred by nature. 
Identifying the true parity partner of a nucleon is also an 
issue. In the mirror models A^(1535) is usually taken to 
be the negative parity state. This choice however fails 
to reproduce the decay width to a nucleon and rj. This 
might indicate another negative parity state lighter than 
the 7V(1535) flO!|, which has not been observed so far. 

In this paper we apply the parity doublet model to 
a hot and dense hadronic matter and study the phase 
structure of a chiral phase transition as well as a liquid- 



gas transition of nuclear matter. The main interest is the 
phase boundary of the chiral symmetry restoration on 
top of the nuclear matter ground state. A light negative 
parity partner will be considered as a phenomenological 
option and we show how the phase boundary lines are 
affected. The phase transitions in large number of colors 
are also discussed. 



2. PARITY DOUBLET MODEL 

The negative-parity nucleon field is introduced such 
that it transforms in the opposite way as the positive 
parity nucleon under SU(2)lx S\J{2)r [1]: 

1p-L->- Rlp-L, 'ip-B,-)- L^p.R, (2.1) 

with L e SU(2)i and R G SU(2)_r, and this allows a chi- 
rally invariant mass term. The Lagrangian of the mirror 
baryons is given by 

C = '(/i+z^?/'+ -I- V^^i^V- +"^-0 ('0-75^'+ - ^'+75^'-) 
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where a, b and gi^ are the coupling constants of the 
mesons to the baryons. The mesonic part follows the 
standard linear sigma model Lagrangian, 
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with the field strength tensor of the omega meson, F^,y = 
d^ijJi, — d^ijj^. The parameters A, /I and e are related to 
the sigma and pion masses and the pion decay constant 



m vacuum as 
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with niTr = 138 McV, /^r = 93 McV and the vacuum 
expectation value of the sigma field ao — fir- The vacuum 
mass of the omega field is m^ = 783 MeV. Other mesonic 
parameters, rria and 174, as well as the chiral invariant 
mass mo will be fixed to reproduce the known properties 
of nuclear matter. 

The mass term of (|2.2p is diagonalized by the mass 
eigenstates of the parity doubled nucleons, N^ and N^ 
via 
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with sinh(5 = ~{a + b)a/2mo. The masses of N± are thus 
given by 
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When the chiral symmetry is restored (ctq = 0), the two 
states become degenerate, 7TiAr+ — rriN- = mo- In con- 
trast, in the naive transformation under chiral symmetry, 
i.e. V'i-L — >• Lip±L and ip±B. ~^ R'4'±R, the nucleon mass 
is generated solely by the spontaneous chiral symmetry 
breaking and thus vanishes with its restoration. 

Applying a mean field approximation '^^j the thermo- 
dynamic potential is 
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where 7^ denotes the spin-isospin degeneracy factor of 
a nucleon and Uf (nj) is the Fermi-Dirac distribution 
function for a fermion (an anti-fermion) defined by 



(2. 



with the energy Ei = ^jp + mf and the effective chem- 
ical potential jj,* = ^ — g^i^o- Imposing the stationary 
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*^ Rotational invariance requires (tJ) 
component (tJo) does not vanish. 



and only the time- 





set A 


set B 


TUN- [MeV] 


1500 


1200 


mo [MeV] 


790 


790 


m^ [MeV] 


370.63 


318.56 


ffw 


6.79 


6.08 


a 


13.00 


9.16 


b 


6.97 


6.35 


fi [MeV] 


199.26 


147.50 


A 


6.82 


4.75 



TABLE I: Set of parameters with gi 
nucleon mass is mM+ = 935 MeV. 



hdj. The vacuum 



condition, dVt/da = dfl/dujo — 0, one obtains the cou- 
pled field equations: 
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The net baryon number density is given by 



p = E T» 



(27r): 



[uf ~ n/] 



(2.10) 



The model parameters are fixed to reproduce the nu- 
clear matter saturation properties at T = 0: 



£'/A(po) - ™Af^ 
po = 0.16fm"^ 



-16 MeV, 



(2.11) 



Originally, the negative parity state has been chosen to be 
iV(1535) [a, 0I- However this choice is incompatible with 
the width of decay into r\N . Also, the mass difference 
between Ai"(1535) and AL|_ is about 560 MeV, that is well 
above the QCD scale Aqcd ~ 200 MeV. The spontaneous 
chiral symmetry breaking might not play much role on 
the physics of such excited states. Another candidate, 
which is lighter and has a stronger coupling to the pion, 
is considered to be a phenomenological option [l3, [ill . 
We will thus consider two sets of model parameters as 
given in [l^j (see Table U): With these sets of parameters 
the incompressibility defined by 
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is 430-510 MeV. This is slightly higher but still compati- 
ble to the suggested range, 200-400 MeV, from the anal- 
ysis of giant mono-pole resonances [I4I . Note that finite 
(74 decreases the incompressibility, e.g. K = 440 MeV 
(tojv- = 1500 MeV)and K = 370 MeV {rriN- = 1200 
MeV) for 34 — 3.8 [lOj- This however does not change 
much either the (pseudo)critical density or the order of 
phase transitions. For simplicity of our calculations, we 
will therefore take 54 = 0. 
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FIG. 1: The ratio of the qq condensate in medium to its 
vacuum value. The calculation was done at T = with the 
parameter set A. For set B there is no significant difference 
from the above result. 



0.1 



CD 

A 

c 

V 



0.01 




0.6 0.8 



1.2 1.4 1.6 1. 



FIG. 2: The expectation value of the sigma field as a function 
of baryon chemical potential at T = 50 MeV. 



3. THERMODYNAMICS OF THE MODEL 

The density dependence of the quark condensate (qq) 
can be extracted using the Feynman-Hellmann theorem 
and it is model independent to the leading order in nu- 
cleon density |13j . 
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with the nucleon sigma term E^v = 45 ± 8 MeV [14 1. 
Fig. [T] shows a comparison of the low-energy theorem 
()3.ip with the in-medium quark condensate calculated 
in the parity doublet model. The model describes the 
correct linearity at low densities up to around the sat- 
uration density po- The value of the condensate, how- 
ever, decreases somewhat stronger than Eq. p. II) . This 
can be considered as an artifact of our simple model 
and taking other parameters (with non-vanishing g^) or 
more systematic treatment would make a better agree- 
ment with the empirical value of Sat. In-medium chi- 
ral condensate has been explored in chiral perturbation 
theory constrained by the nuclear matter properties and 
shown to have a strong pion mass dependence, and a 
substantial deviation from the linear density approxima- 
tion was found already at p ~ 0.2 fm~^ where the ratio 
is ^ 0.7 [l^- In that framework the two-pion exchange 
correlations with virtual A(1232) resonances are respon- 
sible for stabilizing the quark condensate, which are not 
incorporated in the current mean field model. 

It is remarkable that the parity doublet model at T = 
describes a liquid-gas transition of nuclear matter and 
a first-order (set A) or crossover type (set B) of chiral 
symmetry restoration jlO| . At an intermediate temper- 
ature, above these phase transitions, the order parame- 
ter typically shows a double-step structure as shown in 
Fig. [5J The first drop in (cr) appears much below the 



chiral crossover and this is traced back to the remnant of 
the liquid-gas phase transition. The second drop in {a) 
happens at a higher ^b associated with the chiral symme- 
try restoration. The two-step structure is observed only 
at rather low temperatures, and at high temperatures 
{a) eventually experiences a single reduction due to the 
chiral crossover. Note, as mentioned above, the first re- 
duction will weaken when improvements are included to 
reproduce Eq. p.l|) . 

In the parity doublet model one can deal with in- 
medium mesons and baryons explicitly. Besides the 
baryon masses given in Eq. (j2.6p . we define the effective 
meson masses in matter by 

where f7 is given in Eq. (j2.7p . The total meson number 
density is given by 



PM = ^ li 
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where 7^ is the corresponding mesonic degeneracy factors, 
and rib is the Bose-Einstein distribution function. 
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with the corresponding particle energy Ei = ^Jfp- -\- Mf. 
For simplicity, here we consider only a, tt and uj mesons, 
although other mesons, especially p with jp — 9 may 
somewhat contribute too. The meson density depends 
on temperature and it is instructive to compare it with 
the baryon density. We consider the ratio of the meson 
to the baryon number density, 
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to determine a trajectory in the T-fj,B plane. This 
trajectory separates a meson-dominated region from a 
baryon-dominated one and characterizes a meson-to- 
baryon "transition" in dense matter. This is similar 
to the definition adopted in Ref. Il6| where the statis- 
tical model parameters T and fiB were extracted from 
the hadron yield ratios. While the statistical model suc- 
cessfully describes the particle abundances at chemical 
freeze-out, it cannot handle the chiral phase transition. 
In the parity doublet model, on the other hand, the nu- 
clear matter ground state is built in and the chiral sym- 
metry restoration is also described. Thus, the ratio p.Sp 
can be used to find the critical points on the phase dia- 
gram within a self-consistent chiral approach. 

In Fig. [3] we show the phase diagram of the parity dou- 
blet model for the two parameter sets discussed above. 
At zero temperature the system experiences a first-order 
liquid-gas transition at fiB = 923 MeV and the baryon 
density shows a jump from zero to a finite value p ^ Q. 
The order of chiral phase transition and its location de- 
pend on the set of parameters, especially on mass of the 
negative parity state [l3|. Roughly speaking this phase 
transition occurs when the baryon chemical potential 
reaches the mass of the negative parity state, ^b '^ mj^^. 
If we take the most frequently used value tun- = 1500 
MeV, then in addition to the nuclear liquid-gas phase 
transition we obtain a weak first-order chiral transition 
at p ~ 10 po- With a lower mass ttin- = 1200 MeV we 
get no true chiral phase transition but only a crossover 
at much lower density p ~ 3 po ^^ . The liquid-gas tran- 
sition survives up to T = 27 MeV. Above this temper- 
ature there is no sharp phase transition but the order 
parameter is still attracted by the critical point, corre- 
sponding to the first drop in (cr). This makes an addi- 
tional crossover line terminating at the liquid-gas critical 
point. Another crossover line corresponding to the chiral 
symmetry restoration follows the steepest descent of the 
second reduction in (a). With increasing temperature 
the two crossover lines become closer and finally merge. 
The former crossover line, associated with a trace of the 
liquid-gas phase transition, is not well determined from 
the susceptibility when two bumps are getting closer to 
each other. 

In contrast, the trajectory defined in Eq. (|3.5p is ba- 
sically driven by the density effect with the hadron 
masses being not far from their vacuum values. The 
line is almost independent of the parameter set and goes 
rather close to the liquid-gas transition line. The chi- 
ral crossover and the meson-baryon transition lines in- 
tersect at iT,fiB) ~ (150,450) MeV. The parity dou- 
blet model thus describes 3 domains: a chirally broken 
phase with mesons thermodynamically dominating, an- 



other chirally broken phase where baryons are more dom- 
inant and the chirally restored phase, which can be iden- 
tified with quarkyonic matter [ij, [l^- It is worthy to 
note that this point is fairly close to the estimated triple 
point at which hadronic matter, quarkyonic matter and 
quark-gluon plasma may coexist [16j . 



4. CHIRAL CROSSOVER IN NUCLEONIC NJL 
MODEL 

There exist other chiral models which have attempted 
to describe nuclear matter ground state 0, |j]- The 
vector-type interactions at finite density, isoscalar vec- 
tor and a mixed scalar- vector ones, play a significant 
role to reproduce correctly the nuclear saturation point. 
Fermions in those works were introduced in the naive as- 
signment and the nucleon thus becomes massless in the 
chiral limit when chiral symmetry is restored. In this sec- 
tion we give a brief discussion about the thermodynam- 
ics at finite temperatures of a nucleonic Nambu-Jona- 
Lasinio model (NJL) [4] and compare to the result in the 
parity doublet model. 

At pb = the vector interaction does not affect the 
gap equation. 
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(4.1) 

where m% denotes explicit symmetry breaking mass and 
Gs the strength of scalar four-fermion interaction. An- 
other parameter is a momentum cutoff A which regu- 
larizes the integral in Eq. (|4.1|) . The parameters fixed 
to reproduce Eqs. (j2.11|) and (|2.12|) at zero temperature 
are given to be A = 0.4 GeV, Gs = 1.677 GeVfm^ and 
m% = 41.3 MeV [1]. At T = in addition to the nuclear 
liquid-gas transition this model also predicts the chiral 
crossover- type transition at p ~ 3po- Here we check the 
prediction of this model at finite T. Assuming a second- 
order chiral transition in the chiral limit, rriN — >■ for 
771^ — 0, the gap equation determines the critical tem- 
perature as 
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*^ Finite 34 = 3.8, which gives a slightly better fit for the incom- 
pressibility K, does not change its order, but makes a first-order 
chiral transition somewhat stronger |lC| . 



where L2[z] — J2^=i z^/n^ is Euler's dilogarithm. Sub- 
stituting the parameter values one finds for the critical 
temperature Tc = 450 MeV. Numerical calculation with 
explicit breaking leads to a crossover with slightly higher 
temperature Tc = 510 MeV. Obviously this is too high 
to be accepted as a reliable critical temperature compat- 
ible with the Lattice QCD predictions [l^l- On the other 
hand, this Tc is above the cutoff and therefore the model 
can not be used near the transition region. Although 
larger A can cure this drawback, another problem ap- 
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FIG. 3: The phase diagram for the set A (left) and for the set B (right). The crossovers correspond to the local/global maxima 
in the order parameter susceptibility. 



pears, i.e. the bare nucleon mass becomes lower than 
three times the bare quark mass [4] . 

In contrast, the parity doublet model leads to Tc ^ 190 
MeV at which the nucleon and its negative-parity part- 
ner remain massive mN+ ~ rriN- '^^ rno- Therefore, this 
model satisfies the relevant phenomenological constraints 
not only on nuclear matter properties but also on the chi- 
ral crossover at small chemical potential in a reasonable 
way. 



5. PHASE STRUCTURE FOR LARGE N^ 

An interesting issue to address is how the phase struc- 
ture is modified by increasing the number of colors Nc- 
In large Nc limit baryons arise as a soliton in the Skyrme 
model picture. In the current hadronic model we deal 
with the baryons as elementary particles, and thus there 
is no direct access to its dynamical description in the 
large Nc limit, where microscopic structure of the baryon 
is particularly important. Besides, it still remains uncer- 
tain at which scale a Skyrmion collapses (see e.g. Ref. f2(j| 
and references therein). Nevertheless, it is worthwhile to 
study a tendency for various Nc, which may partly bridge 
a gap between the real and large Nc QCD. 

Following Ref. [2l| we adopt the Nc dependences in 
the model parameters as 

mN± -^ (Nc/i) mM± , Too -> (^c/3) toq , 
a -^ ^/Nc/3 a , b ^ y'Njs b , 

9u -> VK/^ 90., u^ y^vs U , 

A^(3/iV,)A, (5.1) 

and the meson masses and the parameter p. do not carry 
Nc dependences. Inserting these scaling factors into the 
gap equations, one can perform thermodynamical calcu- 
lations. In the following, without loss of generality, we 
take zero temperature and consider the liquid-gas tran- 
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FIG. 4: The (pseudo)critical chemical potentials for chiral 
(circle) and liquid-gas (triangle) for various Nc- The negative 
parity state was taken to be 7V_(1200). 



sition and chiral symmetry restoration in cold nuclear 
matter. 

As seen in section |3l how far the two transitions are 
separated essentially relies on the mass of the nega- 
tive parity state: The liquid-gas transition from vac- 
uum to nuclear matter ground state happens roughly at 
fiB ^ 'mN+- When its parity partner is taken to be 
A^(1535), the initial mass gap between the parity part- 
ners, which is generated by the spontaneous chiral sym- 
metry breaking, is rather large ~ 560 MeV. Compared 
with the case taking 7V(1200), more matter compression 
is required in order to decrease the gap and consequently 
to restore the chiral symmetry. This mass gap between 
the parity partners scales linearly in Nc and therefore no 
coincidence of the liquid-gas and chiral phase transitions 
is generally expected. In fact, the actual calculation in 
this model follows this consideration, as shown in Fig. 2] 
With increasing Nc the two phase transitions are system- 



atically shifted toward higher chemical potentials. One 
sees the linear-A'c behavior in hb for around Nc > 10. 
On the other hand, if the effective coupling constants 
increase less slowly with Nc as discussed in Ref. |2(J], 
the two phase transitions may coincide in the large Nc 
limit 22] . More microscopic view of the baryons at large 
Nc might change this prediction, but this study is beyond 
the scope of this paper. 

Deconfinement dynamics which is not dealt with in the 
current model might also modify the above result, espe- 
cially the phase boundaries at higher temperatures. In 
fact, a pure NJL model for constituent quarks predicts a 
chiral phase boundary of elliptic shape for any Nc- While, 
NJL models with Polyakov loops (PNJL) [23 describes a 
hierarchy in quark and gluon sectors with respect to Nc so 
that the deconfinement phase transition always happens 
at the critical value in a pure gluon system, Tq ^ 280 
MeV [13. 



6. CONCLUSIONS 

We have investigated thermodynamics of nuclear mat- 
ter in a parity doublet model within the mean field ap- 
proximation. The model describes the nuclear matter 
ground state at zero temperature and a chiral crossover 
at zero chemical potential at a reasonable temperature, 
which are the minimal requirements to describe the QCD 
thermodynamics. The first-order phase transitions ap- 
pear only at low temperatures, below T ^ 30 MeV. Nev- 
ertheless, at higher temperature they still affect the order 
parameter which exhibits a substantial decrease near the 
liquid-gas and chiral transitions. If the chiral symme- 
try restoration is of first order, criticality around the end 
points of the two first-order phase transitions will be the 
same due to the identical universality class j24| . 

The trace of the nuclear liquid-gas transition is mostly 
overlapped with the meson-to-baryon "transition" char- 
acterizing the ratio of the particle number densities, 
Pm/pb = 1- Its trajectory has only a weak parameter- 
dependence and the intersection to the chiral crossover 



appears at (T,/zb) ~ (150,450) MeV. The chiral phase 
boundary at lower temperature, on the other hand, is 
strongly affected by the negative parity state. The heav- 
ier parity partner induces a larger baryon-rich domain 
in Nambu-Goldstone (NG) phase and thus a strong en- 
hancement in the baryon density is not necessarily ac- 
companied by the chiral symmetry restoration. The true 
parity partner of the nucleon needs to be identified. This 
requires to solve the model dynamically, since in NG 
phase the chiral symmetry does not dictate any non- 
trivial relations between the masses of parity partners 
and their couplings [25]. 

When increasing Nc the two phase transitions are 
shifted to higher chemical potentials and remain sepa- 
rated. This tendency is similar to the observation in the 
PNJL model \V^: The threshold at which the baryon 
number density becomes finite and a first-order chiral 
phase transition do not coincide for any Nc but are close 
with each other. This could indicate that the quarky- 
onic matter [l8| is realized in the chirally restored phase 
with confinement. However, this seems to be unlikely, 
without introducing any gapless modes in a Wigner-Weyl 
phase, since chiral anomalies in low-energy theories do 
not match with those in QCD [2^|. Further investiga- 
tions along with the trace anomaly, which is also respon- 
sible for appearance of a scale in QCD, may bridge the 
thermodynamics of the toy model used in this paper to 
the origin of mg. 
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